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Nonlinearity in macroscopic mechanical system plays a crucial role in a wide 
variety of applications, including signal transduction and processing syn¬ 
chronization and building logical devices However, it is difficult to 
generate nonlinearity due to the fact that macroscopic mechanical systems 
follow the Hooke’s law and response linearly to external force, unless strong 
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drive is used. Here we propose and experimentally realize a record-high 


nonlinear response in macroscopic mechanical system hy exploring the an- 
harmonicity in deforming a single chemical bond. We then demonstrate the 
tunahility of nonlinear response hy precisely controlling the chemical bond¬ 
ing interaction, and realize a cubic elastic constant of 2 x 10^* N/m^, many 
orders of magnitude larger in strength than reported previously^’^^"^®. This 
enables us to observe vibrational bistate transitions of the resonator driven 
by the weak Brownian thermal noise at 6 K. This method can be flexibly ap¬ 
plied to a variety of mechanical systems to improve nonlinear responses, and 
can be used, with further improvements, to explore macroscopic quantum 
mechanics 

Nonlinearity dynamics of a macroscopic mechanical resonator can be mod¬ 
elled by a Duffing oscillator. Its equation of motion under external drive can be 
expressed as: 

mx + jx + kx + ax^ = FdriveCos(a;f). (1) 

Here m, 7 and k are the mass, dissipation rate and linear spring constant of the 
resonator, and ax^ is the Duffing nonlinearity with a the Duffing constant (cubic 
elastic constant). Under weak drive, the nonlinear response is negligible due to 
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its cubic dependence of amplitude x of the resonator, and the resonator behaves 
like a simple harmonie oseillator. This is the well-known Hooke’s law of elastie- 
ity On the other hand, under strong enough drive, nontrivial dynamies of the 
resonator emerges. A famous example is the existenee of Euler instability where 
bifureation oeeurs when the drive strength reaehes a eertain threshold F^. Sinee 
the intrinsie Duffing eonstant in maeroseopie meehanieal systems is weak and 
as Fc is reversely proportional to -y/a, very strong drive is needed in order to ex¬ 
plore the Duffing nonlinearity for applieations This brings in various unfavor¬ 
able faetors: it leads to high power eonsumption and moreover, strong drive 
brings in strong parametrie noise far beyond the Brownian thermal noise even at 
room temperature A method to enhanee substantially the Duffing eonstant 
a for a maeroseopie systems beyond its intrinsie limit is therefore of paramount 
import anee. 

Our system eonsists of a maeroseopie meehanieal resonator tightened to an 
anehor via a ehemieal bonding strueture shown sohematieally in Fig. la. Strong 
nonlinearity is aehieved when the resonator moves in the x-direetion, deforming 
the ehemieal bond. This is beeause, although the resonator alone follows the elas- 
tieity theory with linear dynamies, the response of the ehemieal bond is highly 
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anharmonic. This is confirmed by the chemical bond’s energy curve 17chem(a^) 
(Fig. lb) obtained by density functional theory calculations (see Methods for de¬ 
tails). For simplicity, we have set for the minimum of Uchem{x) (Fig. lb) as a; = 0. 
By applying on the chemical bond an external control force F, the resonator’s 
equilibrium position can be shifted, and the spring constant will be modified by 
Ak = -d^Uchem/dx^. As shown in Fig. Ic, when the resonator’s equilibrium po¬ 
sition is far away from the atom contact, the chemical bonding interaction is weak. 
When we tune the control force, i.e., shift the resonator’s equilibrium position to¬ 
wards a; = 0, the strength of Ak first reaches a local maximum and then drops to 
zero at the place where the chemical bonding attraction reaches maximum. The 
Duffing constant from chemical bonding, i.e., a = changes 

differently (Fig. Id). Its strength reaches a local maximum, then drops to zero be¬ 
fore approaching the maximum attraction point, and finally changes its sign and 
increases in strength. At the maximum attraction point, a = 2 x 10^^ N/m^, while 
A/c is zero. 

Fig. le shows the threshold driving force Fc as a function of the macroscop- 
icity of the resonator (here characterized by the resonator’s mass) with/without 
using the chemical bonding force induced nonlinearity. By introducing chem- 


4 



ical bonding force at the maximum attraction point, can be reduced signifi¬ 
cantly. The term proportional to in the expansion of chemical bonding potential 
Uchemix) will generate a d.c. force as well as a second harmonic term of the drive 
signal, but they contribute nothing to bistate dynamics of our current interests. 
While higher-order terms beyond the Duffing nonlinearity are negligible. 

We demonstrate the above idea by using a macroscopic doubly clamped 
beam whose fundamental vibrational mode couples to a gold-atom contact, as 
shown in Fig. 2a. The atom contact is made by electric current migrations (see 
Methods) on a nano-bridge which anchors the beam to a stiff electrode (Fig. 2b). 
The dimensions of the beam is / x tc x t = 50 fxm x 1.5 /rm x 0.51 /rm. The 
mass of the beam is approximately 0.2 ng. The vibrations of the beam deform 
the gold-gold bonds of the contact along x direction. The beam has an intrinsic 
frequency a;o/27r = 1.58 MFIz, a quality factor Q = 3100, and an intrinsic spring 
constant /cq = 10 N/m. The device is placed in an ultrahigh vacuum chamber and 
has an environment temperature 6 K. 

Fig. 3a plots the measured resonator’s frequency shift Aa;/27r as well as the 
change of spring constant A A; by tuning the control force F. The tuning is done by 
the Lorentz force generated by passing a current through the beam. The resonance 
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frequency of the beam is measured under weak drive, where the response of the 


beam is in linear regime. From the change of spring constant, we estimate the 
Duffing constant as 


1 



/d^Ak 

V dF^ 


k^ + { 




( 2 ) 


with ^ ^ 0.83 the shape constant of beam. In order to reliably obtain the Duffing 
constant from the data with noise, we have smoothed the a using a running av¬ 
erage of 5 points. Fig. 3b plots a as a function of the change of spring constant 
Ak. 


Fig. 3b plots the measured Duffing constant a as a function of Ak. The 
maximum strength of the Duffing constant achieved in our experiments is (2.1 ± 
0.3) X 10^® N/m^, with an enhancement of 4 orders compared to the intrinsic 
nonlinearity due to the elongation of the beam (ao = 2.2 x 10^^ N/m^), where the 
strength of Ak is still very small, i.e., only a few percentage of the intrinsic spring 
constant. When we further shift the beam closer to the electrode, the beam and 
electrode jump to contact, resulting disappearance of the vibration mode. Such in¬ 
stability can be avoided, for example, by exploring a resonator with larger spring 
constant or by lowering the environment temperature (see Supplementary Infor¬ 
mation). In a practical system, apart from the short range chemical bonding force. 
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the Duffing constant can also have contributions from van der Waals (vdW) force 
of the longer range, whose strength is however, small as shown in Fig. 3b. 

We then set the control force at point c in Fig. 3b, and increase the driving 
strength to observe the nonlinear response of the beam. The observed hysteresis 
response is plotted in Fig. 3c. It is different from the hysteresis response of in¬ 
trinsic nonlinearity due to the elongation of the beam, as shown in Fig. 3d. The 
observed hysteresis response in Fig. 3c shows an opposite direction, and the cor¬ 
responding driving force is much weaker (about two orders smaller), indicating 
a negative and much stronger Duffing constant of the beam, as predicted by the 
theory 

In the following we present a demonstration of a dynamical effect of the 
strong nonlinear response. It is known that noise can induce bistate transitions in a 
driven Duffing oscillator However, such transitions in macroscopic mechanical 
resonators are only observed by introducing strong artificial noise that is much 
large than the Brownian thermal noise at room temperature, due to weak intrinsic 
nonlinear response We are able to observe in our system the bistate transitions 
activated by the Brownian thermal noise even at cryogenic temperatures. In doing 
so, we tune the control force to a point with the Duffing constant approximately 
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— lx 10^^ N/m^ and drive the system to bistable state regime. By reeording the 
amplitude of vibrations of the beam we observed a switehing behavior, as shown 
by the green trajeetory in Fig 4a. 

To verify that sueh switehings are indeed from the Brownian thermal noise 
of the beam, we introduee an amplitude modulation to the driving signal with am¬ 
plitude as Fdrive(f) = -^drive + 5Fdrive COs(f2f) ,with = 0.5 Hz and (JFdrive = 
0.18 pN. The amplitude of vibrations of the beam shows periodie switehings 
(Fig. 4a, purple trajeetory) instead of random switehings. Fig. 3b shown the eorre- 
sponding power density, S'mod(^) with modulations and S'noise(^) without modu¬ 
lations, from whieh we define a signal to noise ratio: SNR(H) = S'mod(fi)/*S'noise(^) 
Then, by using standard stoehastie resonanee theory (see Methods), we estimate 
the total foree noise as = (3.6 ± 0.6) x 10“^® N/\/Hz. It agrees nieely 

with the resonator’s Brownian thermal noise, estimated by = AmuJoksT/Q 
whieh equals to 3.3 x 10“^® N/-\/Hz. In the limit H approaehing zero, the mod¬ 
ulation beeomes a perturbative foree signal ()FdriveCos(a;f) added to the driving 
foree with the same phase. In the rotating frame of the driving signal, the bistate 
dynamies is deseribed by an over-damped double well whose shape is tuned by 
(iFdrive and beeomes sensitive near the bifureation point. By applying a weak 
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force ^-Fdrive = 2 fN that is only several times larger than the resonator’s Brownian 
thermal noise foree, signifieantly ehange in the statisties distribution is observed, 
shown in Fig. 4e. The sensitivity of sueh response to external foree is limited by 
the total foree noise \/deviee, besides the thermal noise, paramet- 
rie noise from the eireuit also eontributes to the total foree noise, with the weak 
power density being We estimate that is 11 dB below 

the level, eorresponds to a parametrie noise temperature of 500 mK, and is 
mainly limited by our room temperature deteetion eireuit (see Methods). 

In eonelusion, we have demonstrated a highly eontrollable nonlinearity in 
a maeroseopie meehanieal system, with its fiuetuation dynamies dominated by 
Brownian thermal noise. The universal existenee and the small seale of the ehem- 
ieal bonding foree make our method applieable to eurrent widely-used miero- and 
nano-meehanieal systems in improving their nonlinear responses. With further 
improvements, nonlinearity indueed quantum behaviors in maeroseopie meehan- 
ies are foreseeable in the type of systems deseribed here. 
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Methods 

Theoretical description of the system. Density functional theory calculations 
were employed to estimate qualitatively the chemical bonding force. Due to sta¬ 
bility considerations, a small structure of two gold-atom clusters was adopted in 
the calculation (see Sec. 1 of Supplementary Information). By changing only 
the relative position x of the left cluster relative to the right one while otherwise 
keeping the relative positions of all the gold atoms fixed, we obtained a chemi¬ 
cal bonding energy function Uchera{x), from which we extracted the nonlinearity 
by calculating the fourth order derivative of the total energy. Due to the exist¬ 
ing systematic error in performing density functional theory calculations, the sec¬ 
ond order derivative of the total energy shows slightly an oscillatory behavior at 
large distance range. We smoothed the calculated data by fitting it to the func- 
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tion In practical case, however, the relaxation of gold atoms ean oeeur 

whieh will ehange the energy - displaeement x eurve. We have eonsidered sueh 
effeets and found out that the results do not ehange signifieantly (see Sec. 1 of 
Supplementary Information for detailed data and comparison). To estimate the 
long-range forces, we considered analytically the van der Waals (vdW) attraetion 
and the eleetrostatic force for the geometric model of two gold cylinders close 
to eaeh other, similar to the ease of our nano-bridge strueture. To estimate the in- 
trinsie maeroseopie meehanieal property, we modeled the maeroseopie resonator 
as a doubly elamped beam. The lowest vibrational mode has a weak intrinsic 
nonlinearity due to elongation. In calculating the threshold drive Foree in Fig 
le, we took the meehanieal quality faetor 3000 and the ratio between the beam’s 
thiekness and length t/l > 0.005. 


Fabrication of the sample. The deviee was fabrieated on a eommercial Silieon- 
On-Insulator wafer using nanolithography. The nano-bridge conneeting the beam 
and the stiff eleetrode is about 200 nm in length and 80 nm in diameter. Foeus- 
ion beam (FIB) was adopted to narrow it down to less than 50 nm. Once sueh 
a deviee, i.e., the doubly elamped beam with a suspended nano-bridge, was sue- 
eessfully fabrieated, we plaeed it in the ultrahigh vaeuum environment of the eryo- 
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genic system, and then electro-migrate the bridge to form an atomic point contact 
After the atom contact was produced, the beam was disconnected from the 
stiff electrode and its lowest vibrational mode was measured. The position of the 
equilibrium of the beam is controlled by applying a d.c. current through the beam 
so the atom interaction can be tuned. The separation of the atom contact can be 
controlled with a precision about 1 pm without any feedback control for thou¬ 
sands of seconds (see Sec. 2.4 of Supplementary Information for the data). We 
applied a voltage bias less than 200 mV between the beam and the stiff electrode 
to compensate the contact potential of the atom contact so that the contribution of 
electrostatic force in our experiments can be negligible (see Sec. 3 of Supplemen¬ 
tary Information). We decreased the separation of the atom contact gradually by 
increasing the control force, while monitoring the frequency of the resonator. In 
doing this, we make the beam vibration amplitude around 50 pm until a jump to 
contact happens, and after that, a finite resistance (typically around 1 Kf2) can be 
detected between atom contact which shows good contact without contamination. 
Then we decrease the control force until the beam detaches from the stiff elec¬ 
trode. The typical force required to detach the beam is several nN, indicating that 
the atom contact is a small gold-atom cluster 
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Device characterization and electrical parametric noise. The Brownian ther¬ 
mal noise on the beam is 3.3 x 10“^® N/a/Hz, and the thermal motion amplitude 
noise on resonanee is about 1.0 x 10“^^ m/A/Hz aoeordingly. The intrinsic non¬ 
linearity is measured from the frequency response to the driving strength in a 
standard way based on the Duffing nonlinear model 

To do this, we have pulled the beam far apart from the stiff electrode with 
contact separation larger than 30 nm, so the interaction due to the atom contact 
is negligible. The intrinsic Duffing constant obtained is 2.2 x 10^"^ N/m^ (see 
Sec. 2.5 of Supplementary Information for data). To estimate the electrical para¬ 
metric noise in our system, we model the mechanical resonator as LCR elements 
and analyze the noise in the electric circuit (see Sec. 3 in Supplementary Infor¬ 
mation). The electrical noise in the circuit produces a current noise /noise passing 
through the beam, which results in a force noise Bllaoise of power density 

In our experiment, there are mainly three sources of electrical noise. The 
first one is the Johnson-Nyquist noise S'^ = ARksT, with T the room temper¬ 
ature, which produces a force noise with power density The second one is 
the current leakage from the input of the voltage preamplifier We model the 
preamplifier back-action by a current source similar to that described in reference 
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In doing this, we assume that there is no eorrelation between voltage impreei- 
sion and the baek-aetion so 5^^ leads to foree noise with power density 
Another one is from the r.f. driving signal, whieh generates a phase noise 
and works as a equivalent foree noise on the beam with power density S'p^a- We 
estimate the total parametrie foree noise as: and in the 

experiments, = 0.9 x 10“^® N/a/Hz, with power density about 11 dB 

below the 5^^ (see See. 3.2 of Supplementary Information for data). 


Measurement of total force noise on the beam. The dynamies of our system are 
modeled by the standard Duffing oseillator of equation of motion 

mx + 'yx + kx + ax^ = FdrWeit) COs{ut) + Fnoise(t), (3) 

where the dissipation rate 7 = muo/Q, the spring eonstant k = mu^, Tdrive(^) = 
-Fdrive + ^Tdrive cos(f2f) is the amplitude modulated driving foree, and Fnoise(^) 
is the total noise with power density The minimum foree that drives the 

system into the nonlinear regime where bifureation oeeur is Fc, with the eorre- 
sponding vibration frequeney being 2ttuJc, and the amplitude Xc ean be ealeulated 
from Eq. (3) Near the nonlinear bifureation point, we transform this equa¬ 
tion to an over-damped one by following the standard procedure For the case 
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of modulation frequency f2 being much smaller than the decay rate ojqIQ of the 


system, the signal to noise ratio (SNR) is related to the system’s noise power as 



(4) 


with Xm the half of the vibration amplitude of the bistable states, Su the nonlin¬ 
earity frequency induced shift from linear resonance peak, and 7 ^ the measured 
random switching rate without modulation. So from the measured SNR we ob¬ 
tained 
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Figure 1: Concept of the system and theoretical results, a, A macroscopic 
resonator tighten to an anchor via chemical bonds. The displaced resonator can 
compress (up) and stretch (down) the chemical bond of gold-atom contact. b,c,d, 
Density functional theory calculation results of (b) the chemical bonding inter¬ 
action energy Uchem{x), (c) the modified spring constant Ak = —d'^Uchem/dx'^, 
and (d) the enhanced nonlinearity coefficient a = — {l/Q)d'^Uchera/dx'^ as func¬ 
tions of the resonator displacement x. e, Estimated threshold drive force Fc as a 
function of the resonator mass, m, with: (i) intrinsic nonlinear response of the res¬ 
onator (dark blue), (ii) enhanced nonlinear response of the resonator (light blue) 
by chemical bonding interaction, as indicated in gray line in (b,c,d). where chem¬ 
ical bonding induced linear response Ak = 0. (See Methods for details of the 
model.) 
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Figure 2: Experimental setup, a, Scanning electron microscopy of the sample 
in false color. The macroscopic resonator is a doubly clamped silicon beam with 
thin layer of gold deposited on it, with dimension I x w x t = 50 fxm xl.5 fxm 
xO.51 /im and total mass approximately 0.2 ng. The center of the beam has 
horizontal displacement (x). In the presence of a 6 T external magnetic field along 
the z direction, the electric current (/) can excite and detect the motion of the 
beam, with the schematic circuits shown, b. Nano-bridge connecting the beam to 
a stiff electrode, before experiments, c. Cartoon plot of the atom contact generated 
on the nano-bridge indicated by ‘c’ in b, the gold-gold bonding interaction is then 
tuned by force F which is controlled by a d.c. current through the beam, and the 


electrostatic interaction of the contact is minimized applying a d.c. bias on the tip. 
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Figure 3: Tuning the macroscopic nonlinear response by chemical bonding 

force, a, Frequency shift (left axis) and the corresponding change in effective 

spring constant Ak (right axis) as a function of the external force F that pushes 

the atom contact close to the beam, to do this, frequency resonance under weak 

drive with the beam in near linear regime are measured (inset), b. The tunable 

chemical bonding interaction can simultaneously change the Duffing nonlinear 

response a and the effective spring constant hy Ak (minus is for convenience). 

The measured Duffing constant (blue dots) is attributed to the short range chemical 

bonds (orange line), the long range van der Waals interaction (blue thick line), and 

the intrinsic Duffing constant due to the beam elongation (gray dashed line), c, d, 
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Hysteresis response of the system measured at the point indicated by ’c’ in b, and 
measured with beam pulled far apart (> 30 nm) from the stiff electrode so that 
only intrinsic nonlinearity due to the beam elongation contributes (d). 
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Figure 4: Brownian thermal noise induced bistate transitions, a, Switching be¬ 
tween bistable states with constant-amplitude drive (upper panel) and modulated- 
amplitude drive (lower panel). To observe bistate transitions, the driving fre¬ 
quency is fixed in the middle point of hysteresis response and the drive amplitude 
is Fdrive = 4 pN, and the modulation frequency and the amplitude are = 0.5 
Hz and ^Fdnve = 0.18 pN, respectively, b, Power spectrum density. Compared to 
the constant-amplitude drive (upper panel), the amplitude-modulated drive (lower 
panel) induces an additional peak in the power spectrum density at the modulation 
frequency = 0.5 Hz. c. The amplitude distribution of the bistable resonator de¬ 
pends on the driving amplitude. Very different histograms for constant drive with 
.Fdrive = 4.000 pN (Left panel) and 4.002 pN (Right panel) are observed. 

The change in amplitude distribution can be understood using an effective double 
well potential (sketched in the inset). The total electrical parametric noise is 11 dB 
below the Brownian thermal noise in our experiment. 


























